Strong coupling phases of the t — t' Hubbard model 
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We study the strong coupling regime of the t — t Hubbard 
model, filled up to the level of the van Hove singularities, by 
means of an exact diagonalization approach. We characterize 
the different phases of the model by the different sectors of the 
Hilbert space with given quantum numbers. By looking for 
the ground state of the system, we find essentially the com- 
petition between a state with incipient ferromagnetism and 
other mimicking a d-wave condensate, which has the lowest 
energy in a large region of the phase space. 
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During recent years there have been increasingly ac- 
curate measures by angle resolved photoemission spec- 
troscopy of copper-oxide compounds, giving much in- 
sight into the phenomenology of these materials!. Near 
the optimal doping for superconductivity the hole-doped 
compounds use to show extended van Hove singularities 
close to the Fermi level, which, are located near the high- 
symmetry points (n, 0), (0, 7r)BB. On the other hand, in 
the carrier free regime the materials show antiferromag- 
netic correlations, with a dispersion relation which has 
peaks at the points (±7r/2, ±7t/2)H. A most interesting 
problem is therefore to understand the drastic change 
that the Fermi surface may suffer by the influence of 
dopinga. 

The framework that has been proposed to address such 
theoretical issues is that of the t — t' — U jpodeld (or its 
strong coupling version, the t—t 1 — J modeE), as it is gen- 
erally believed that strong correlation effects have to be 
responsible for the electronic properties of the cuprates. 
Next-to-nearest-neighbor hopping t' has to be introduced 
for a more accurate description of the dispersion relation 
in the insulating phasal The distinctive feature of the 
t — t' — U model is that the level of half- filling does not co- 
incide with the level corresponding to the two van Hove 
singularities. Thus, it should be possible to establish 
a clearer separation between the effects of the antiferro- 
magnetic correlations and those due to the appearance of 
the extended saddle-points. There have been attempts 
to propose a purely electronic mechanism of supercon- 
ductivity in systeps with van Hove singularities close to 
the Fermi levelEflL3. What is essential in those models 
is the existence of some enhanced channel favoring the 
exchange of singlet pairs. They represent an alternative 



to the picture earlier proposed in which the pairing in- 
teraction is supposed to apse from the short-range anti- 
ferromagnetic correlationstil. 

In the present paper we study the correlations which 
may dominate at the van Hove singularities, by perform- 
ing the exact diagonalization of the t — t' Hubbard model 
in a 4 x 4 lattice (with periodic boundary conditions). 
The hamiltonian of the model is 



H=-tJ2 



ct„c 



+ u 2l n ^ n a (!) 



where a =|, j, the first sum is over nearest neighbors i, j, 
the second sum over next-to-nearest neighbors, and n, is 
the electron number operator at site i. A typical contour 
map of the dispersion relation (for t' < 0.5 ) is shown 
in Fig. 1. We are specially interested in the situation in 
which the van Hove shell, comprising the four degenerate 
states at (tt, 0) and (0, it), is half-filled. We will pursue 
the determination of the lowest energy state in each of the 
sectors with different quantum numbers. By looking for 
the ground state of the model we will be able to study the 
interplay among the different phases of the system and, 
as long as we are relying on a quantity (the ground state 
energy) that probes the lattice as a whole, we may hope 
to predict properties with less influence from finite-size 
effects. This prpcedure is similar to that applied to the 
Hubbard mode£3, where it appears to be safer than the 
evaluation of correlation functions on the small lattice 
scaleM 




FIG. 1. Contour map of dispersion relation for the t — t' 
Hubbard model. The van Hove shell (VHS) is shown. 
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Our starting point is the lattice in Fig. 1, filled with 
12 particles for the 32 available one-particle states. For 
the noninteracting theory (U — 0), this means that the 
ground state is sixfold degenerate. The six states are 
given by the different occupations of the van Hove points 
A at (tt, 0) and B at (0, tt). We may classify them accord- 
ing to the total momentum P and total spin S, ending 
up with a S = 1 triplet at P = (tt, tt) 

|r 1 )= c + T c+ T |0) 

\To) = ^ (4 T 4i + 4i4 T )l°) 

|T_ 1 )= C + i 4jO) (2) 
a spin singlet at P = (tt, tt) 

\P) = ^(4Al- c M^)\0) (3) 
and two spin singlets at P — (0, 0) 

I^> = ^(4t4|-4 t 4i) \o) 

|s> = ^(4t4i + 4 t 4i)|o) (4) 

The states \D) and \S) (as well as \T ) and \P)) differ in 
the quantum number associated to the transformation by 
exchange of the two components of the momentum, that 
is one of the generators of the point symmetry group. 
The triplet state would be the precursor of a state with 
nonzero magnetization for systems with larger size, while 
\P) would correspond to a paramagnetic or antiferromag- 
nctic state with no tendency to the uniform alignment of 
spins. On the other hand, the \D) state mimics the fluc- 
tuation of pairs between the two van Hove singularities, 
providing a simplified version of a <i-wave condensate. In 
what follows we study how the degeneracy between these 
states is broken when the interaction is turned on, inter- 
preting the corresponding ground state as an incipient 
signal of what should be the dominant correlation in the 
model. 

By considering the interaction of the particles in the 
van Hove shell with the rest of closed shells in the Fermi 
sea (but without allowing yet for particle-hole screening 
processes) we see that the above ground state degeneracy 
is partially lifted. To first order in U, the energy of each 
of the above multiplets is 

E lT) =E kin + ^(^ + 2)^f + 0(U>) 

E\ D) =E km + l((^f + l) 2 -l) + 0(^) 

E lP) =E km + ^((^f + iy + ?j+0(Ui) 

%> = E km + ^ ((^f + I)' + l) + 0(U 2 ) (5) 



where E^m is the kinetic energy of the particles, tlfs is 
the number of particles in closed shells and N is the num- 
ber of lattice sites. We see that to this perturbative order 
the states |T) and \D) have both the lowest energy. If we 
had an extended saddle point, for instance, and the pos- 
sibility of placing a number to of particles near A or B in 
a more degenerate level, it is clear that aligning the spin 
of the m particles would produce a state with an energy 
~ E kln + (U/N) (n FS /2 + to) n FS /2 + 0(U 2 ), while any 
state made of a condensate of singlet pairs could hardly 
lower its energy from w E kin + (U/N) ((nps + m ) /2) 2 + 
0(U 2 ). It seems therefore that, in the weak coupling 
regime of a model with a high density of states close 
to the Fermi level, a state with a macroscopic amount 
of magnetization is energetically favoredtll One has to 
bear in mind, however, that this argument works only 
for weak coupling constant. The influence of particle- 
hole processes and screening turns out to be crucial in 
the presence of the van Hove singularities. In our 4x4 
lattice model one may check that, already at values of 
the coupling constant as small as U « 2t, there is a value 
of t' below which the state corresponding to \P) dressed 
with particle-hole processes gets lower energy than that 
of |T). Our exact diagonalization study becomes rele- 
vant mainly in the strong coupling regime, where there 
is no forecast of what may be the symmetry of the ground 
state. 

As an illustration of the effects at large U we have plot- 
ted in Fig. 2 the evolution of the minimum values of the 
energy in the different sectors with the respective quan- 
tum numbers of the \T), \D) and \P) states, in the region 
of small t' and U — 8t. From the computational point of 
view, we have reduced the Hilbcrt space in each case by 
using all the generators of the point symmetry groupEj, 
and we have implemented a Davidson algorithm to ob- 
tain the lowest energy state in each sector. The plot in 
Fig. 2 shows that above a certain value of t' the ground 
state is found in the sector with the symmetry of the 
\D) state and, in fact, it corresponds to the original \D) 
state in (^) dressed with a large amount of particle-hole 
processes from the Fermi sea of our model. However, in 
the region of very small t' (< 0.02i), the lowest energy 
is reached in the sector of the \P) state. This is consis- 
tent with what we expect from the Hubbard model, that 
is what we recover in the limit t' — > although with a 
significant deviation from half-filling. The ground state 
that we find in that regime corresponds to what should 
be a state with paramagnetic or antiferromagnetic prop- 
erties in larger lattices. It is worthwhile to remark that 
such state becomes a highly excited state in our study 
for sufficiently large values of the next-to-nearest neigh- 
bor hopping, being also overriden by the lowest energy 
state with spin S — 1 above some value of t' . The two 
crossings of levels that we appreciate in Fig. 2 have a 
genuine quantum character, since we have already seen 
that at very weak coupling constant the ground state has 
to become essentially either the triplet state in @ or the 
\D) state in (|). 
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FIG. 2. Minimum value of energy in the sectors of |ZJ),|T) 
and \P) states (asterisks,diamonds and triangles resp.) in the 
region of small t' and U = 8t. 

We have also studied the behavior of the model at in- 
creasing values of U and next-to-nearest neighbor hop- 
ping t' = 0.3i, that is adequate for making contact with 
the phenomenology of the cuprates. The sectors which 
may have candidates for the ground state of the system 
have again the quantum numbers of \T), \D) and \P), at 
least up to U = 12t. It can be checked that states with 
momentum different from (0, 0) or (tt, tt), or spin different 
from or 1, have much higher energy. Furthermore, any 
state with the quantum numbers of IS 1 ) is always an ex- 
cited state of the system. In the mentioned range of the 
coupling constant, the ground state turns out to be the 
\D) state conveniently dressed by particle- hole processes 
from the closed-shells of the Fermi sea. We have repre- 
sented in Fig. 3 the difference in energy of such state 
with the state of minimum energy in the sector of |P), 
and in Fig. 4 with respect to the corresponding lowest en- 
ergy state built from \T). It becomes clear that the state 
which could signal the appearance of antiferromagnetic 
order has no chance of being the ground state of the sys- 
tem, for the value of t' = 0.3t considered. Moreover, the 
results plotted in Fig. 4 show that the state which mim- 
ics the cf-wave condensate has always lower energy than 
the triplet state up to the largest value U = lOt consid- 
ered. It is important to note that at the largest values 
of U the energy difference between the two levels has a 
very slow decline, what makes unlikely that the triplet 
state may become the ground state for a near value of 
U . The situation seems similar to the case of Nagaoka 
ferromagnetism, where this effect is found in the small 
4x4 lattice for very large values of U ( > 60< ) , although 
the nature of the ferromagnetism may be quite different 
in the present model. The persistence of the \D) state 
as the ground state of the system over a large region 
of the parameter space provides good reason to believe 
that the corresponding d-wave condensate may continue 
playing that role in the larger lattices. 
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FIG. 3. Difference between minimum energies of the sec- 
tors containing \D) and \P) states (t'=0.3). 
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FIG. 4. Difference between minimum energies of the sec- 
tors containing \D) and \T) states (t'=0.3). 

The possible crossing of levels between the singlet with 
d-wave symmetry and the S = 1 state is also reminiscent 
of a phenomenon studied in small clusters of the hexag- 
onal lattice. In those systems the regime in which the 
singlet prevails as the ground state has been linked to 
the effect of pair binding above half-filling. The hexag- 
onal lattice at half-filling has strictly two Fermi points, 
what makes likely that such pair binding may rely on the 
resonance of electron pairs at the Fermi level. The im- 
portant difference with respect to the model with the two 
van Hove singularities is that the density of states for the 
hexagonal lattice goes to zero at the Fermi points, which 
is just the opposite of the situation that we are facing. 
This may explain why the ground state with pairs fluc- 
tuating between the two van Hove points is more stable 
in this case, persisting over a much wider range of values 
of U than in the clusters studied in Ref. 

The results that we have obtained have to be taken 
anyhow with the reserve inherent to the use of a small 
lattice in the exact diagonalization approach. As we have 
pointed out before, our analysis is likely to be relevant in 
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the study of the strong coupling regime, where the effects 
of screening due to particle-hole processes are difficult to 
assess by any other method. Moreover, our results ap- 
ply directly to the zero temperature regime of the sys- 
tem, while the use of other approaches like the quantum 
Monte Carlo method make very difficult to perform ex- 
trapolations to the limit of vanishing temperature. This 
drawback seems to be overcome anyhow in the study of 
the t — t' Hubbard model of Ref. [l?], where it is claimed 
that at U = At and t' = 0A7t a signal of ferromagnetism 
is present in the model. For those values of the cou- 
plings our picture does not differ much qualitatively from 
what is shown in Figs. 3 and 4, but this is consistent 
with the fact that the signal found in Ref. |l7| is so weak 
that requires lattices as large as 16 x 16 for the ferro- 
magnetism to be observed. It has to be stressed that 
the value t' = QA7t is very close to the point in which 
the Fermi sea degenerates into a pair of straight lines, 
making the analysis of the scattering processes of Ref. 
|l7| in terms of ladder diagrams quite appropriate. Our 
approach, however, has to be relevant for a completely 
different regime, that arises when particle-hole processes 
induce strong screening effects in the modeltS. 

Of course, in an ideal analytical framework one should 
be able to deal simultaneously with the description of 
screening processes and scattering in the particle-particle 
channel. One way to accomplish this would be to adopt 
a renormalization group approach, taking first into ac- 
count the scaling of the different channels and study- 
ing then the behavior of the response functions. The 
two channels that are not irrelevant, apart from the for- 
ward scattering channels, are the scattering V of pairs 
about the same van Hove singularity and the scatter- 
ing V of pairs from one singularity to the other. The 
regime in which the intra-singularity screening dominates 
making V < V (what happens above certain value of t') 
seems to lead to an instability favoring the condensation 
of pairs with opposite amplitude in the two van Hove 
singularities".. This is nothing but a new version of the 
Kohn-Luttinger mechanismc3, in which the anisotropy of 
the screening effects arise in a natural way. The super- 
conducting instability may therefore be competing with 
the above mentioned ferromagnetic instability, at least 
for values of t' not very close to 0.5t. It would be worth- 
while to develop some computational scheme with which 
to perform an extrapolation of our results to larger lat- 
tice sizes, in order to see whether the evidence found for 
d-wave superconductivity survives in the strong coupling 
regime. 

We want to thank F. Guinea, R. Hlubina and S. Sorella 
for useful discussions during the development of this 
work. 
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